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ABSTRACT 

This model assumes the baryons, radiation, three famihes of massless neutrinos, and cold dark matter 
were mutually thermalized before the baryon number was fixed, primeval curvature fluctuations were 
subdominant, and homogeneity was broken by scale-invariant fluctuations in a new dark matter compo- 
nent that behaves like a relativistic ideal fluid. The fluid behavior could follow if this new component 
were a single scalar field that interacts only with gravity and with itself by a pure quartic potential. 
The initial energy distribution could follow if this component were gravitationally produced by inflation. 
The power spectra of the present distributions of mass and radiation in this model are not inconsistent 
with the measurements but are sufficiently different from the adiabatic cold dark matter model to allow 
a sharp test in the near future. 

Subject headings: cosmology: dark matter — galaxies: formation — cosmic microwave background 



L INTRODUCTION 

One motivation for the search for alternative models for 
structure formation is that we arrived at the commonly 
discussed adiabatic cold dark matter (aCDM) picture after 
just a few false starts. This might be because the early 
universe is simple enough that there are only a few ways 
structure could have originated, or because we were lucky, 
or perhaps because more than one model is viable at the 
present level of constraints. It seems prudent to continue 
the search for viable alternatives before we learn whether 
they are needed. 

The candidate presented here draws elements from a 
phenomenological model (Hu 1999) that allows an accept- 
able fit to the measured power spectra of distributions of 
matter and radiation, and a model with a physical prove- 
nance within the inflation scenario (Peebles 1999a) but a 
poorer fit to the measurements. Our new picture has the 
observational advantage of the former and a physical basis 
that simplifies the latter. It has elements in common with 
the phenomenological models analyzed by Bucher, Mood- 
ley, & Turok (1999), but crucial differences that make the 
present model viable. 

We start with the idea that, since the dark matter may 
interact only weakly with ordinary matter and radiation, 
some or all of it may interact only with gravity (Peebles & 
Vilenkin 1999a; 1999b, hereafter PVb). Such dark matter 
would be gravitationally produced, as a squeezed state, by 
inflation (Ford 1987; Grishchuk & Sidorov 1990). 

We discuss initial conditions from inflation in the next 
section, evolution of the departures from homogeneity in 
§3, and tests of the power spectra of the matter and radi- 
ation in §4. 

2. INITIAL CONDITIONS 

The dynamical components are the cosmic microwave 
background (the CMB), three families of massless neutri- 
nos, baryons, cold dark matter (CDM), and a new dark 
component that acts like an ideal fluid with the equation 
of state Py = Py/2i. The primeval energy density contrasts 



satisfy 



5~t = 5u = -Sh 
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The last part expresses the isocurvature condition, where 
p is the total mass density in which the baryons and CDM 
are initially subdominant. The power spectrum of 5y is 
nearly scale-invariant: k^Pyik) is initially close to con- 
stant. The first part of equation (|^) says the fiuctua- 
tions in the usual matter components are adiabatic. This 
can follow if the chemical potentials of the neutrinos and 
CDM vanish and all these components are in mutual ther- 
mal equilibrium that is broken after the baryon number 
is frozen and before the CDM is nonrelativistic. The rel- 
ativistic fluid behavior of py can follow from a field that 
interacts only with gravity and with itself by a quartic 
potential, with the action 



S= I a^d^xdtiy,,y'y2-\yy4), 



(2) 



and energy density py — y'^/2 + (Vy)^/2 + Xy^/'i. 

When the frequency of the field oscillation is large com- 
pared to the Hubble parameter a/ a, equation (||) expressed 
in conformal time t = J dt/a{t) is the action in Minkowski 
coordinates for y = ay. Since the energy of y in Minkowski 
spacetime is conserved, the mean energy density in y scales 
as Py oc a{t)~'^ (Ford 1987, Peebles 1999b). This means py 
can be large enough to serve as a primeval seed for struc- 
ture formation but remain small enough not to interfere 
with the standard models for light element production and 
gravitational structure formation. 

Fluctuations in py are well approximated as linear acous- 
tic waves from the end of infiation, when the field starts 
oscillating, through the characteristic acoustic oscillation 
time divided by the density contrast (Peebles 1999b). The 
acoustic wave model fails with the appearance of features 
that resemble shock waves. If the scale- invariant spectrum 
of y extends to small wavelengths these shock-like features 
appear well before the field fluctuations of interest to as- 
tronomy appear at the Hubble length. The analysis in Pee- 



bles (1999b) indicates that this does not spoil the acoustic 
wave model on larger scales. 

If the field y in equation (||) exists it will have been 
excited, with a near scale-invariant spectrum, by inflation 
(Ford 1987). Kofman & Linde (1987) considered the near 
classical evolution of y in inflation when the potential for 
the inflaton cj) also is quartic, V{(j)) = X^(j)'^/4. We assume 
this same eternal inflation model. 

We assume the dimensionless parameter in equation (pi) 
is in the range 

A0 < A < 0.01. (3) 

The lower bound makes the y-field energy subdominant to 
the inflaton during inflation (PVb); otherwise y assumes 
the role of the inflaton (Felder, Kofman & Linde 1999). 
The upper bound is chosen so y is close to constant across 
the present Hubble length. This follows from a consider- 
ation of the competition between the freezing of quantum 
fluctuations that tend to drive the field value away from 
zero and classical dissipation as the field rolls to the min- 
imum of its potential at y = 0. Early in inflation these 
processes are close to statistical equilibrium. As the value 
of the Hubble parameter H decreases equilibrium eventu- 
ally is broken, at expansion parameter a^, and y thereafter 
evolves almost as a classical field. Under the upper bound 
in equation (0) the expansion from Og to the end of in- 
flation is large enough that y is close to constant across 
our Hubble length. [| We assume that at our position the 
field value at a — a^ is close to the characteristic value at 
equilibrium, A (?/'*) = 3H'^/8it^ (Starobinsky & Yokoyama 
1994). The perturbations to y added from a = ag to the 
end of inflation produce a near Gaussian scale-invariant 
fluctuation spectrum with variance per logarithmic inter- 
val of wavcnumber (PVb eq. [39]) 



SyY ^_ A 

y J 6tt^ logg ax/ Up 



(4) 



As discussed in the next section the large-scale pertur- 
bation to the present thermal background radiation (the 
CMB) is 5T/T = {Sy/y)/5f, where / is the ratio of mass 
densities in radiation and neutrinos to the energy in the 
j/-field. The fit to the observed temperature variance per 
logarithmic interval of I, (ST/T)'^ = (1 x 10~^)^, requires 



A« 1 X lO^V • 



(5) 



The value of / depends on the model for the origin of 
ordinary matter and radiation; we must treat / as an 
adjustable parameter. The standard model for the light 
elements requires / ^ 10, leaving a small window of con- 
sistency between equations (^) and (||). 

3. EVOLUTION 

The evolution of the distributions of matter and radi- 
ation in linear perturbation theory is computed by the 
usual methods. We discuss only aspects that differ from 

^This follows by adding numerical factors to equation (33) in 
PVb, to get the y-field relaxation time Ty = iTy/EH~^X~^'^, with 
Ty = H^ at a = Oe . This with the expansion factor Oa, /ap ~ e^" to 
the end of inflation from the time of freezing of the fluctuations we 
see, with the condition a^ < Qp, fixes the bound on A. 



the usual casejj 

The dynamics of the fluctuations are governed by two 
events; the transition from radiation- to matter-dominated 
expansion, and Hubble crossing, when relativistic stress 
gradients and gravity have comparable dynamical effects. 
We consider first a large-scale mode that crosses the Hub- 
ble length after matter-radiation equality. The residual 
entropy fluctuation, 



<7 = 5c- -(5rcl 



5c = -5^ 



(6) 



where p^e\ — P-y + Pi^ + Py, becomes important near matter- 
radiation equality and before pressure causes the mode to 
oscillate. On these large scales the relative distribution 
of the y-field and the familiar radiation species is irrele- 
vant because their gravitational effects exactly cancel be- 
fore Hubble crossing and are negligible afterward. The 
relation between the present distributions of the cold dark 
matter and the CMB shows an interesting effect. The 
gravitational potential in the matter-dominated regime is 
related to the initial entropy fluctuation as $ = a/5, and 
$ is related to the density perturbation by the Poisson 
equation (Kodama & Sasaki 1986). The CMB anisotropy 
due to the gravitational redshift is 
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This is the same as the adiabatic CDM (aCDM) model, 
and different from the isocurvature CDM (iCDM) model in 
Peebles (1999a) where primeval fluctuations in the CDM 
are balanced by the CMB and neutrinos and (5T/T|grav = 
— 2$. It is helpful to the construction of a viable model 
to assume the initial photon distribution follows that of 
the species that is responsible for gravitational structure 
formation in the matter-dominated epoch (Hu 1999). This 
puts the photons that initially are hottest where the gravi- 
tational potential becomes the deepest, so the temperature 
fluctuation is reduced as the photons move out of the po- 
tential well. The consequence is that in aCDM and the 
present model the observed ratio of matter to radiation 
fluctuations follows from a near scale-invariant primeval 
fluctuation spectrum, while in the iCDM model a fit to 
this ratio requires a substantial tilt to increase small-scale 
fluctuations over large. 

The evolution of small-scale fluctuations in the CDM 
that cross the Hubble length before matter-radiation equal- 
ity is surprisingly sensitive to the behavior of the new 
y-component. At Hubble crossing, stress gradients in 
the relativistic components cause them to oscillate. The 
CMB density oscillates as an acoustic (sound) wave with 
amplitude given by the initial conditions in such a way 
that observationally acceptable peaks result from scale- 
invariant initial conditions. The y-component in our 
model also oscillates as an acoustic wave. Aside from 
the neutrinos, this keeps the radiation distributions al- 
most balanced. The CDM amplitude appears at the Hub- 
ble length after moderate growth from the initial value 
(Fig. |l|a). As in aCDM, the modest further increase of 

■^Though there is no gauge ambiguity in our initial conditions 
numerical stability in the evolution requires careful choice of gauge 
and metric variables (Hu 1999). 
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Fig. 1. — Time evolution of Fourier amplitudes in a mode with 
very short comoving wavelength for (a) our perfect fluid model for 
the y-component and (b) a free massless particle model. Note 
the zero crossing of the CDM amplitude (5c in (b). In a mode at 
k ~ 0.15/1 Mpc" ^, (5c crosses zero near matter-radiation equality, 
producing the zero of P(k) in the dashed line in Fig. 2. 



5c to matter-radiation equality leaves the usual k~'^ sup- 
pression of small-scale power and an observationally ac- 
ceptable present mass fluctuation spectrum from initially 
scale-invariant fluctuations. If instead the ^/-component 
were a free gas of massless particles, free streaming would 
cause an imbalance with the acoustic oscillation of the 
CMB, as in Figure Ep. The resulting metric perturbation 
reverses the sign of Oc- The reversal can only be produced 
during the radiation-dominated epoch, so there is a mode 
near the Hubble scale at matter-radiation equality that 
is caught in the act of reversal, producing a zero in the 
linear power spectrum, and spoiling the fit to large-scale 
structure measurements. 

4. PHENOMENOLOGY 

Figure compares our model predictions to all signif- 
icant measurements of the CMB temperature anisotropy 
and to the Peacock & Dodds (1994) compilation of mea- 
surements of the power spectrum of the galaxy distri- 
bution. We assume a scale-invariant initial spectrum of 
fluctuations in y (as in eq. 0), standard recombination, 
J7m = 0.35, r^A = 0.65, h — 0.8, and two values of the 
baryon density, fifc/i^ = 0.02 and fife/i^ ^ o.012. The for- 
mer is close to the central value of Buries et al. (1999) 
based on the deuterium abundance D/H = 3.4 x 10^^. 
Kirkman et al. (1999) consider the most secure bound to 
be D/H < 6.7 x 10"^; this abundance scales the baryon 
density to the lower number. 

Following Tegmark (1999) and Miller et al. (1999), we 
find a crude estimate of x^ for the CMB temperature 
anisotropy by treating all data points as independent with 
Gaussian distributions of errors. The first two lines of Ta- 
ble 1 list the reduced x^ employing all the data (58 points), 
the selection in Miller et al. (1999) ("A"; 24 points), and 
the remaining data plus the COBE DMR resuhs ("B"; 42 
points). The third line is the best fit aCDM model from 
Tegmark (1999). Values for the full data set and selec- 
tion "B" are arguably less secure because they are based 
on a more heterogeneous set of methods. The ~ 10% 
calibration uncertainty, which is not included in these x^ 
estimates, is a serious general barrier to the interpreta- 
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Fig. 2. — Measured and modeled power spectra of the CMB 
and galaxy distributions. The shading of the CMB error boxes is 
according to area and corresponds to the la errors X the width 
of the experimental window; cross hatching denotes measurements 
included in set "A". Our models are plotted as solid curves; the 
dashed curve represents a model where the seed field is a gas of free 
massless particles rather than a fluid. 



tion in terms of formal measures of significance. Within 
the calibration uncertainty our low baryon density model 
seems viable, although challenged by the D/H measure- 
ments (Kirkman et al. 1999), while our high density model 
is challenged but we believe not ruled out by the CMB 
measurements . 



TABLE 1 

Approximate CMB x^/^^ 



Model 


AU 


A 


B 


Qbh' = 0.012 
Qbh^ = 0.02 
aCDM 


2.6 

2.7 
2.5 


1.5 
2.0 
1.2 


1.3 
1.3 

1.4 



We get satisfactory agreement with the power spectrum 
of the galaxy distribution. The normalization implies a^ = 
0.84 - 0.86, consistent with the bounds 0.74 ^ as ^ 1.1 
implied by the abundance of rich clusters of galaxies at 
our model parameters (e.g. Viana & Liddle 1999). 

Our model differs from those of Bucher et al. (1999) 
in two ways. First, the CDM density perturbations fol- 
low the CMB (eq. [|l[). We noted that this suppresses 
the large-scale anisotropy of the CMB, allowing a near 
scale-invariant primeval spectrum and making the peaks 
in the anisotropy spectrum at .£ > 100 more prominent. 
The same effect follows from a coherent superposition of 



the CDM-isocurvature and neutrino isocurvature modes of 
Bucher et al. (1999); one cannot use a linear combination 
of the individual power spectra. Second, motivated by 
equation (g), we model the isocurvature departure from 
homogeneity by a component that behaves as a perfect 
fluid rather than a gas of free massless particles. One sees 
in Figure 1 and the dashed curves in Figure 2 that the 
CMB anisotropy is not much affected but the free parti- 
cle model produces a zero in the mass power spectrum (in 
linear perturbation theory) at an undesirable wavelength. 
The power spectra of the present distributions of matter 
and the CMB depend on the cosmological parameters in 
different ways from the aCDM model. The heights of the 
peaks in the CMB spectrum depend on fl^h^ in opposite 
ways (Hu & White 1996): here the odd- numbered peaks 
represent rarefaction of the photon fluid in the potential 
wells and hence decrease when the baryon density is in- 
creased, as one sees in Figure 2. Also, since our model has 
no initial metric fluctuations whose decay in the radiation- 
dominated epoch enhance the peaks, the peak values do 
not increase with decreasing $7^^^- The lesson here is 
a general one: cosmological parameters derived from a 
model fit are provisional no matter how securely fixed 
within the model until the model is unambiguously es- 
tabhshed. 

5. DISCUSSION 

We conclude that our model is viable but likely to 
be critically tested by CMB anisotropy measurements in 
progress. The same is true of the aCDM model, of course. 

Our model can be adjusted; here are four considera- 
tions. First, we use isocurvature initial conditions. There 
may be a significant adiabatic perturbation from the in- 
flaton, or, in other inflation models, from the stress of the 
j/-field fluctuations during inflation. Second, we place A in 
a narrow window (eqs. |3| and (5)). If A ^ 0.01 the fluctu- 
ations in py at the end of inflation have positive skewness, 
so the primeval fluctuations in the CDM mass distribu- 
tion are non-Gaussian with negative skewness. Models 
with positive skewness are seriously constrained (Frieman 
& Gaztafiaga 1999); negative skewness may be interesting 
for structure formation. The second moments needed for 
the tests in Section 4 have not been analyzed for this case. 
Third, one may ask whether some or all of the CDM is 
in fields that interact only with gravity and themselves 
by potentials that would have to be more complicated 
than the quartic considered here. PVb present prelimi- 
nary elements of a model for this more complicated case. 
Fourth, we have assumed standard recombination. One 
could imagine stars present at z ~ 1000 delay the rapid 
drop in ionization; that would shift the peaks in the CMB 
spectra to smaller £ and substantially change the signifi- 
cance of this test. 

The field y is a new hypothesis. Its parameter A is not 
exceedingly small, however, and by moving the seed for 
structure formation from the inflaton we remove the re- 
quirement for a specific value of the very small parameter 
A^. But closer consideration of these issues might best 
await observational developments. 

It may be significant that the structure formation his- 
tory in our model is only mildly different from aCDM. 
Perhaps this is telling us viable phenomenological models 
already are limited: they have to approximate aCDM. Or 



perhaps our imagination in exploring concepts like gravi- 
tationally produced matter is limited. 
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